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Abstract 

In this note we consider the inviscid limit for the 3D Boussinesq equations without diffusion, 
under slip boundary conditions of Navier's type. We first study more closely the Navier-Stokes 
equations, to better understand the problem. The role of the initial data is also emphasized in 
connection with the vanishing viscosity limit. 
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1 Introduction 

The aim of this note is to study the L°°(0, T; L 2 (Q)) convergence, as v vanishes, of (Leray-Hopf) 
weak solutions to the 3D Navier-Stokes equations (NSE), towards smooth solutions of the 3D Euler 
equations. The dependence of the rate of convergence in terms of different hypotheses on the initial 
data is studied and the results are also applied to handle the problem of convergence of solution of 
the 3D Boussinesq equations to those of the Euler-Boussinesq. 

We first study the problem with constant density (set for simplicity equal to one) and then, 
in the final section, we treat the Boussinesq equations. In particular, we start by considering the 
inviscid limit for the NSE in a bounded domain Q C H 3 with smooth boundary F := <9f2 ^ 0. 

For the reader's convenience we recall that when T is non empty, for the NSE with Dirichlet 
boundary conditions (and with v > 0) 

d t u u - uAu" + (u v ■ V) u v + Vp u = infix (0, T], 

V-u u = inOx(0,T], 

u u = onTx(0,T], 

11^(0, x) = Uq in Q, 
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in general one cannot have convergence (even in weak norms) towards smooth solutions of the Euler 
equations, even with the same initial data Uq = u E 



(1) 



d t u E + (u E ■ V) u E + Vp E = in Q x (0, T], 

X7-u E = inOx(0,T], 

u E -n = onrx(0,T], 

u E (0,x)=u E inU, 



see e.g. the review in Constantin [9] and Mazzucato |21| . In fact, even if both u E and u v are 
very smooth and both exist in [0, T] (for some positive T independent of the viscosity) certain 
extra- assumptions are needed in order to show, at least, that 

as v -> 0+ u v (t) -> u E (t) in L 2 (0), uniformly in t G [0, T). 

Some necessary and sufficient conditions, related with the dissipation of energy in a boundary-strip 
of width depending on v, have been detected by Kato j!3j . See also recent developments in Temam 
and Wang [25 1, Wang |26| . and Kelliher [13]. The lack of convergence is due to the boundary layer 
created from the difference between the tangential velocity of the Navier-Stokes solution and that 
of the Euler solution at the boundary: The first vanishes, while we do not have control on the 
tangential velocity of the Euler equations. 

Better results can be obtained in the case of the NSE with Navier's boundary conditions. In 
Iftimie and Planas [12] it is considered the following initial-boundary value problem 



d t u v - vAu u + (u u ■ V) u u + X7p u = in Q x (0, T], 

V-u u = inOx(0,T], 
(2) u"-n = onrx(0,T], 

[V(u v )n + Pu v ] t!m = onrx(0,T], 
u u (0, x) = Uq in O, 

where T>{u u ) = \ [Vu" + (Vu^) T ] is the deformation tensor, (3 > is a constant (the friction co- 
efficient) and \D(u v )n + /3n^] ta n is the tangential component of the vector T>(u u )n + (5u v . This 
system is very close to that originally proposed by Navier [22] and studied analytically (in the sta- 
tionary case) starting from Solonnikov and Scadilov |23] , In particular, the Navier's slip conditions 
read as {vT>(u u )n + /3u u ]t an = 0, hence in |12] the authors are implicitly assuming the Maxwell 
scaling [20j . with the friction parameter /3 depending linearly on the viscosity. More details on the 
role of Navier's boundary conditions especially for numerical simulations, and some of the crucial 
differences between the two dimensional and three dimensional case, can be found in the review 
paper [7j. In the 2D setting the problem is slightly less-hard and classical results employing slip 
boundary conditions are those of Yudovich [32], J.-L. Lions [16], and Bardos [I]. Interesting results 
in the 2D case are those in [181 [T9] . 

A recent vanishing viscosity result in the 3D for system ([2]) is the following one, (see Theorem 1 
in [12]). 



Theorem 1.1. Let Q be a bounded smooth open set in R 3 and let uq £ H^(Q) a divergence- free 
vector field tangent to the boundary. For each v > consider u^ € L 2 (f2) a divergence free vector 
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field tangent to the boundary such that Uq — > uq strongly in L 2 (Q), as v — > 0. Let u u be a weak 
solution of the Navier-Stokes equations ([2]) with Navier's boundary conditions and with initial datum 
Uq. Let u E G C([0, T]; H 3 (Q,)) be the unique solution of the Euler equations (pQ), with initial datum 
uq, for some T G]0,T max [, being T max the maximal time of existence of the smooth solution of the 
Euler equations. Then, u v converges to u E strongly in L°°(0, T; L 2 (f2)), as v — >• 0. 

Remark 1.2. By inspecting the proof, one can observe that, if the initial data converge in L 2 (Q) fast 
enough, then the same argument implies that sup ig [ 0i T] ll n ^(^) ~ uE {t)\\ 2 = 0{y) and J Q \\Vu v (t) — 
Vu e (t)\\ 2 dr = 0(1). The result here is independent of the parameter /3 > 0. 

Our aim is to study the convergence under some different slip-without-friction boundary con- 
ditions, involving the vorticity. More precisely we will study the following initial-boundary value 
problem 

d t u v - uAu v + (u v • V) u" + Vp" = in Q x (0, T], 

V-u u = inOx(0,T], 

(3) u u -n = onrx(0,r], 

curlu 1 ' x n = on F x (0,T], 

1^(0, x) = Uq in f2. 

and we will show how the convergence-rate can be improved. 

The interest for these vorticity based Navier's boundary conditions is increasing, especially after 
the recent results by Xiao and Xin [291 EH] j Beirao da Veiga and Crispo O HI [SJ |6] and [8] concerning 
strong solutions and strong convergence. See also the related work by Xin et al. [271 EH E] . Next, 
we point out that the connection between the two Navier's type conditions, ([2]) versus ([3]), is 
expressed by the following identity, valid for all tangential vectors r on the boundary T: 

(4) t ■ t = — (curlu x n) ■ r — vu ■ — on T. 
Here t is the Cauchy stress vector defined by 

n 

t(u,p) := n ■ T(u,p) = ^T ik {u,p)n k , 

k=l 

with T(u,p) := —lp + vT>{u). The vector identity (jU) valid on T shows that the two Navier's 
conditions are essentially the same in the case of a domain with flat boundary. Moreover, in a 
general domain they differ by a lower order term. 

Results similar to the present ones have been also recently obtained by Xiao and Xin |30| . 
while the results presented here are part of the Ph.D. thesis of the second author, completed at 
the end of 2011. We point out that the approach in |30| is slightly different, focusing more on 
H 1 ^) convergence in the case of "well-prepared" initial datum Uq = u E G i? 3 (0) for both the 
NSE and Euler equations. In addition, the boundary value problem for the Navier-Stokes equations 
is a non-standard one, slightly different from ([3]). On the other hand, the aim of our result is to 
show how the initial data affect the convergence of the vanishing viscosity even in the energy norm. 
Moreover, we give a simple proof of the the convergence in the case of well-prepared initial data. 
In particular, the main result of this paper is Theorem 13.21 which shows improved convergence in 
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the L 2 -norm and convergence also of first derivatives for the system ([3]), when the initial datum 
has vorticity vanishing at the boundary. Observe that in Theorem 13.11 the vanishing vorticity is 
requested only for the initial datum u E of the Euler equations. The initial data Uq of the NSE 
are just divergence- free vector fields, tangential to the boundary, and converging merely in L 2 (n) 
(hence without any control on the vorticity) to the datum of the Euler equations. 

In the final section we use the results obtained for the NSE equations to tackle the following 
problem: We study the convergence of (v u ,p u ), solution of the viscous Boussinesq equations 

d t v u - i/A/ + («" ■ V) v u + Vq v = -p u e 3 infix (0, T], 

d t p u + (v v ■ V) P v = in Vl x (0, T], 

V-v u = inOx(0,T], 

(5) v v -n = onrx(0,T], 

curlu" x ri = onrx(0,T], 

v u (0,x) = v% inQ, 

pr(o,x) = (% inn, 

toward those of the Euler-Boussinesq equations 

d t v E + (v E ■ V) v E + Vq E = -p E e 3 in O x (0, T] , 

d t p E + {v E ■ V) p E = in U x (0,T], 

V-v E = infix(0,T], 

(6) Tp 

v E -n = onrx(0,T], 
v E (0,x)=v E inn, 
p E (0,x)=p E inn, 

in the energy space L°°(0,T; L 2 (0)), where e 3 is the third vector of the canonical basis in R 3 . 

Plan of the paper In Section [2] we briefly recall the notation, some vector identities, and the 
existence results for the NSE (|3|). Next, in Section [3] we prove two different vanishing viscosity 
results for the NSE, showing the critical dependence on the initial datum. Finally, in Section [3] we 
study the vanishing viscosity limit for the Boussinesq system. 



2 Preliminaries 

We consider a bounded domain n C R 3 with smooth boundary T, say of class C 4 , and n denotes 
the exterior normal unit vector on T. We will use the classical Lebesgue spaces (L 2 (f2), || . ||) and 
(L 2 (r), || . ||p) and the Sobolev spaces (H k (n), \\ -\\k) for k € N (we do not distinguish between 
scalar and vector valued functions). We will denote by (H S (T), \\ . \\ s ,r) the standard trace spaces 
on the boundary T. We will also denote by C generic constants, which may change from line to line, 
but which are independent of the viscosity and of the solution of the equations we are considering. 

We first start by recalling the precise notion of weak solution for the NSE we will use. 
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Definition 2.1. We say that u v G L°°(0, T; L 2 (Q)) n L 2 (0,T; H 1 ^)), weakly divergence-free and 
tangential to the boundary, is a (Leray-Hopf) weak solution of the Navier-Stokes equations ([3]) if 
the two following conditions hold: 



[ [ ( - u v 4> t + vVu u V<j) - {u u ■ V) <(> u v ) dxdr + v [ f u u ■ ( Vnf ■ <f> dSdr 
Jo Jn Jo Jr 

UQ(j)(0) dx, 



(7) 



for all vector-fields <p G C£°([0, T[xf>) such that V^ = 0fflOx[0,T[, and <p ■ n = on T x [0,T[; 
the following energy estimate 



v\\2. 



(8) ~h v (t)f + u fwVvTfdr + v f I u v • (Vn) T ■ u v dSdr <hu ol] 

2 Jo Jo Jr 2 

is satisfied for all t G [0, T] . 

With this definition we have the following result. 

Theorem 2.2. Let be given any positive T > and Uq G L 2 (£l) which is weakly divergence- free 
and such that Uq ■ n = on T. Then, there exists at least a weak solution u u of the Navier-Stokes 
equations ([3|) on [0, T]. 

The proof of global existence of weak solutions in the sense of the Definition 12.11 can be found 
for instance in |29} § 6] . We observe now that our definition of energy inequality is slightly different 
from that in the above reference and we explain now the equivalence. To this end we recall the 
following formulas for integration by parts (see Ref. [2] for the proof). 

Lemma 2.3. Let u and <f> be two smooth enough vector fields, tangential to the boundary T. Then 
it follows 

- I Au <p dx = / VuVcj) dx - (u xn)(pdS + / u- ( Vnf ■ cf) dS, 
Jn Jn Jr Jr 

where uj = curlu. Moreover, i/V- u = 0, then —Au = curl curl u, and 

/ cutIuj (j)dx = — / Aucj)dx= \ w(curl0) dx + / {ojxn)<j)dS. 
Jn Jn Jn Jr 

Constructing weak solutions by the usual Galerkin method as in , we get for the approximate 
solutions {n^jmgKf the following identity: 

- I |u^(t)| 2 dx + v I I | curl-u^J 2 dxdr = — f \uQ m \ 2 dx t & [0,T], 
2 Jn Jo Jn ^ Jn 

where we used as test function the function u 1 ^ itself and the second integration by parts formula 
from Lemma 12.31 with u = (ft = u v m (here the boundary integral vanishes due to the vorticity-based 
Navier's boundary conditions). Then, the usual compactness tools imply that u u m converge as 
m — > +oo to a weak solution u v and the lower-semi-continuity of the norm implies that 



- / \u u (t)\ 2 dx + u / / | curing 2 dxdr < — j \uq\ 2 dx, 
2 Jn Jo Jn 2 J n 
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which is the energy inequality (6.12) in [29]. Finally, by using the first integration by parts formula 
from Lemma 12.31 we get (JSj) - This inequality will be used later on to let some of the calculations 
(which will be otherwise only formal) completely justified. 

To conclude, we recall a well-known existence theorem for smooth solutions of the Euler equa- 
tions (HJ) in Sobolev spaces. 

Theorem 2.4. Let be given u E G H 3 (Q) such that V • u E = and u E ■ n = on F. Then, there 
exists a positive time T = T(\\u E > such that a unique solution of ([I]) exists in 

u E G C{[0,T];H 3 (n)). 

The proof in the case of a bounded domain can be found in Ebin and Marsden [10J and 
Temam [24j . In particular T > „ , for some C > independent of the solution, and in the 

sequel T will be any positive time strictly smaller than the maximal time of existence T max . 



3 Proof of the convergence results 

We start by showing the basic convergence result, which is the counterpart of [X2^ Thm. 1] in our 
setting. 

Theorem 3.1. Let Q be a bounded smooth open set in IR 3 , and let u E £ H 3 (ft), be a divergence- 
free vector-field tangential to the boundary. Let u E G C([0, T]; H 3 (Q,)) be the unique solution of 
the Euler equations (pQ) ; with initial datum uq and defined in some interval [0, T\. Let u v be a 
weak solution of the NSE §3§ with divergence-free and tangential to the boundary initial datum 
Uq G L 2 (ft), and with the vorticity-based Navier's conditions. Suppose also that 



3 , 



K -^11 = 0(1/3). 

Then, su Pie[cvr] \\u"(t) -u E (t)\\ 2 = 0{A) and j Q T \\Vu v {t) -Vu e {t)\\ 2 dr = 0{A). 

Proof. The proof is simply obtained by taking the difference of the equation satisfied by u v with 
that for u E , multiplying by u := u v — u , and integrating by parts over f2 x (0,T). Unfortunately, 
this cannot be done in a so straightforward manner since u u is a weak solution, hence using directly 
u (having the regularity of u u ) as test function is not allowed. We need to pass to an integral 
formulation and to use the energy inequality ([8]) to make the argument rigorous. The reader 
well-acquainted with the argument can go directly to the formula (]14p . 

We first observe that since u E is a smooth solutions of the Euler equations (H|) in [0,T] x Q, 
u E is allowed as test function for the NSE. Then, after certain integration by parts, we get for all 
t G [0, T] 



0) 



f u v {t)u E {t)dx+ [ [ (vVu v Vu E -u v uf + {u v -V)u v u E )dxdT 
Jn Jo Jn > 

+v f f u u {Vn) T u E dSdr = [ u E u u dx. 
Jo Jr Jn 



A further identity is obtained by multiplying the Euler equations by u v . Since u E is a local smooth 
solution everything is well-defined and we get 

(10) / / (uf u u + (u E ■ V) u E u v ) dxdr = f u E u u Q dx. 

Jo Jn Jn 
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Next, by multiplying the Euler equations by u E and by the usual integrations by parts we get the 
energy conservation 



(ii) 



dx = — 
2 



u E \ 2 dx, 



where we used the fact that u E is smooth, tangential to the boundary, and divergence- free. 
Then, by adding together ([8J)- (fTT]) and subtracting (|Uj)- (fTUj) . we get 

rt 



\u(t)\\< 



+ v 



(12) 



o Jn 



\7u u Vu dxdr + 



o Jn 



+ v 



[ [ u v {Vn) T udSdT < 
Jo Jr 



(u ■ V) u E u dxdr 

um 2 



Let us consider the second term from the left-hand side of (|12p . By using the parallelogram equality 
we get 



t 

o Jn 



Vu v VudxdT 



1 



Vu u \ 2 dxdr + 



o Jn 



I Vnt 2 dxdr 



o Jn 

t 



\Vu E \ 2 dxdT. 



o Jn 



Then, we estimate the other two terms from the left-hand side of (|12p . We handle the third term by 
using the higher-order regularity ofu E , namely we use u E G C([0, T];H 3 (Q)) ^ C([0, T]; W 1,00 (Q)), 
to get 

r i r r i 

\u(t)\\ 2 dT, 



[ f (u- V)u E udxdr <C I 
Jo Jn Jo 



for some C independent of v. Concerning the boundary term in (|12p . by using trace theorems and 
Young inequality we get, due to the smoothness of T, 



o Jr 



u v {Vn) T udSdT 



< v I \u u \ |Vn| \u\ dS 



< C ^||n^||r||^||r 

< Ci/||n i/ || 1 /2||'"||i/2 

< Ci/||« l '||5||Vu I '||5|kl|i||V'u||^ 



< -||V«"|M + — NVtilr + CVIKIIIMI- 
4 4 

In particular, to handle the i? 1 / 2 (r)-norm and to remove the zero-order term we have used the fact 
that for functions tangential to the boundary the Poincare inequality holds true (see for instance, 
Kozono and Yanagisawa [15]) and also that from the energy inequality both \\u u \\ and \\u E \\ are 
bounded. 

By collecting all the estimates, from (|12p we get that 



u{t)\\ 2 + u / \\Vu(t)\\ 2 dr < \\u(0)\\ 2 + C 



(13) 



We can now use Gronwall lemma, obtaining 



\u(t) \\ 2 dr + v 



u — u 



.E\\2 



\L°°(0,T;L 2 (n)) 



O(u) 



and 



V^-V^||| 2(0)T;L2(n)) = O(l) 
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exactly as in |12j. This is not the result stated in the Theorem 13.11 but the calculations are given 
to better understand the differences/improvement. 

We make now some slightly different manipulations, in order to show the better rate-of-convergence 
stated in Theorem 13.11 To this end, in (|12p we treat the second term from the left-hand-side as 
follows ^ ^ 

/ / Vu v VudxdT = j I \Vu\ 2 dxdr + / / Vu E Vudxdr. 
Jo Jn Jo Jn Jo Jn 

We arrive now at (cf. (]12p ) 

»t r rt 



(14) 



o 



Vufdr- j j (u-V)u E udxdr + i> [ [ u u (Vn) T udSdr 
Jo Jn Jo Jv 

<-u t [ Vu E VudxdT+ IK0)I ' 2 



'o Jn 2 

We handle the third term from the left-hand-side exactly as before and, by integrating by parts 
the term from the right-hand-side (by using the second identity from Lemma 12. 3p , we get 

-v I I Vu E VudxdT 



o Jn 

v \ \ (u E x n)udSdT + v \ u E ■ (Vn) T • udSdr + v / Au E udxdr. 
Jo Jr Jo Jr Jo Jn 

Then, we get the following inequality 

+ „ /' || V „|| 2(iT < JK21!! _ v f I ul y n f uiSiT 

j (u ■ V) u E udxdr + v / / Au E udxdr 
o Jn Jo Jn 

t 

E 



— v I I (oj x n)u dSdr. 
Jo Jr 

We estimate the absolute value of the space integral from the right-hand-side by using Schwartz 
inequality, trace inequalities, and the regularity of u E as follows: 



(to x n)udS 



< Cv\\u\\^\\Vu\\^ < Cv% +C \\u\\ 2 + -||Vw| 



J u-{Vn) T -udS 



< Cv\\u\\l < Cv\\u\\\\Vu\\ < Cv\\u\\ 2 + -||Vn| 



Au E u dx 



<C[\\ 



u\\ + v 



Then, by using also the energy inequality, we obtain the following differential inequality (cf. 
with ([131)) 



+ v / ||Vu(r)|rdT < ||w(0)|r + C 
Jo 



\u(r)\\ 2 dT + v 2 + 
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By using Gronwall-Lemma we have that 



K " U E \\l°°(0.T;L*(n)) = °^~ 2 ) and H V ^ " Vu E \\ 2 L2( „ T .r 2(n)) = 0(U2), 



ending the proof. □ 

We want now to show better convergence, and this happens if the initial datum belongs to a 
particular sub-class. In particular, we use the same observation made in [8j [30] to show strong 
convergence up to second order derivatives. We prove now the main result of the paper. 

Theorem 3.2. Let Q, be a bounded smooth open set in R 3 ; and let u E £ H 3 (Q,), be a divergence- free 
vector field tangential to the boundary, and such that 

(15) uj e (x) = o Vxer. 

Let u E E C([0, T]; H 3 (Q,)) be the unique solution of the Euler equations (pQ) ; with initial datum u E 
and defined in some interval [0, T\. Let u u be a weak solution of the Navier-Stokes equations ([2]) 
with a divergence free and tangential to the boundary initial datum Uq 6 L 2 {VL) such that 

ho-u E \\ = 0{y). 

Then, 

sup \\u" - u E f Lao(0TL 2 m = 0(v 2 ), and \\Vu u - Vu S ||| 2{0 T i2(n)) = <D(v). 

t£[0,T] 

A critical point in the proof is that of ^having solution to the Euler equations with vanishing 
tangential component of the vorticity, as the Navier-Stokes equations," that is oj e x n = on 
r x [0, T\: In this way one can better estimate the term — f^Au E udx involved in the previous 
calculations. 

In general the boundary conditions for the vorticity cannot be enforced for the Euler equations. 
In addition, also if the initial datum is such that (uj e x n)|p = this is not enough to have the 
same boundary behavior for all positive times. As observed in [3j, by using the vorticity equation, 

wf + (u E -V)u E = (cj e -V)u E , 

by taking the exterior product with the normal unit vector on V, and finally by using that co E x n = 
implies that uj e x n = on T, one obtains that an extra-compatibility condition, generically false, 
should be satisfied by the initial velocity u E . In particular, this implies that the Navier's type 
condition does not persist for positive time and hence excludes the chance of a vanishing-viscosity 
limit in topologies such that the vorticity uj e has Sobolev traces at the boundary. On the other 
hand, by using the fact that for the Euler equations the vorticity is transported by the velocity u E 
and stretched by Vu E , one can employ the well-known representation formula for classical solutions 

(16) oj E {X{a, t),t) = V a X{a, t) co E {a, 0), 
where the path-lines X : Q, X [0, T] — > Q, C R 3 solve the Cauchy problem 

f £ t X(a,t) = u E (X(a,t),t), 
\ X(a,0)=a, 
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for t G [0,T] and for a £ Q. Since u E ■ n = on the boundary, the path-lines starting on the 
boundary remain on the boundary for all positive times. The fundamental effect for our studies is 
the following: Let be given a G T, then X(~a, t) G T for all t G [0,T] and consequently 

uj E (X(a, t),t) xn = [V a X(a, t)u E (a, 0)] x n, V (a, t) G T x [0, T]. 

Being generically the matrix V Q X(a, i) not a multiple of the identity, this implies that generically 
co E (X(a,t),t) xn/0. In particular uj e x n may become non-zero, as soon as V a X induces a 
rotation along any axis not parallel with the normal unit vector passing through a. 

On the other hand, this observation makes clear that it is possible to have another type of 
persistence for the vorticity, by suitably restricting the class of initial data. In fact, the same 
representation formula with the path-lines shows also the following result. 

Lemma 3.3. Let u E G C([0, T]; H 3 (Q)) be the unique solution of the Euler equations If 
oj e {x) = for all x G T, then 

oj E {x,t) = for all couples (x,t) G T x [Q,T]. 

Remark 3.4. The class of solutions with vanishing-vorticity at the boundary we will employ is in 
some sense optimal, since in [5$ it is shown that if u>o • n ^ 0, then there exists a time T > such 
that the oj e xn/0 for < t < T. Moreover this class of initial data in non-empty, since smooth 
and divergence- free functions with compact support satisfy the assumptions. 

Remark 3.5. In the case of the Euler equations with a non-zero external force f the same 
approach shows that a formula similar to (|16j) holds true: 

u E {X(a,t),t) = V a X(a,t)uj E (a,0) + / curl f E {X(a, a), a) da, 

Jo 

hence a sufficient condition to have persistence of the vanishing vorticity at the boundary is ui E = 
curl f E = on T. 

This observation on the persistence of the boundary values for the vorticity points out that 
probably the "w-based" boundary conditions of Navier's-type are much better behaved (in the con- 
text of vanishing-viscosity) than the classical "P(n)-based" Navier's ones. The Navier's conditions 
as in ([2j) involve the symmetric part of the gradient. One can recall that the symmetric part of the 
gradient T>(u E ) := — - - has the following evolution equation 

E£^l + V (u E f + 0\u E ) = -Hn E , 

where as usual & is the derivative along path- lines, 0(u E ) := ) 1 is essentially the vor- 

ticity (since 0(u E )h = |w E x h for each vector h) and Htt e is the Hessian of the pressure. It 
seems that (contrary to the results for uj e or equivalently O) the evolution of the matrix T>(u E ) 
cannot be handled, since the pressure does not disappear and one cannot employ directly a trans- 
port/stretching argument. Hence, the problem related with vanishing- viscosity under the Navier's 
boundary conditions seems to require different tools, even if the friction parameter j3 vanishes. 
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Proof of Theorem \3.2l We can improve a little bit the rate of convergence of Theorem 13.11 by 
assuming that the initial datum is such that the vorticity vanishes at the boundary. In fact, we 
have seen that by assuming (|15p we have co E (x,t) = for all (x,t) £ T x [0, T]. Writing again the 
same energy estimates we employed before to get (fT4l) , and by using the identity 

\7u E \7udxdr = v I I u E ■ (Vn) T -udSdr + v [ [ Au E udxdr, 
Jo Jr Jo Jn 



o Jn 



(notice that now the boundary integral f r (w E X n)u dS vanishes) we can re-do the same calculations 
as before starting from (JT3J) to obtain now 



u{t)\\ 2 + v \\Vu\\ 2 dr < \\u(0)\\ 2 + C / \\u{t)\\ 2 dT + V 



Finally, by GronwalFs inequality we get 

0(v 2 \ 



sup \\u u — u E 1 1 



te[o,T] 

concluding the proof. 



L°°{0,T;L 2 (Q)) 



and HVti" - Vu E \\ 2 L2 ^ T . L 2 m = 0{v), 



□ 



4 Vanishing viscosity for the 3D Boussinesq equations 

In this section we pass to consider the Boussinesq system © and a first step concerns the existence 
of weak solutions. In particular, since the problem is without diffusion in the equation for the 
density, a proper notion of weak solution is needed. Observe that the unknown p u entering the 
equations is not the density, but represents the small variations of density from the constant state, 
hence there is no need to assume p u > 0. 

Definition 4.1. We say that the pair u u G L°°(0, T; L 2 (Sl))nL 2 (0, T; H 1 ^)) and p u <G L°°(0, T; L 2 (Q)) 
is a weak solution of the Boussinesq equations ([5]) if the two following conditions hold: 

[ (-v u (p t - P v ^t + vVv u Vcj) - (v w ■ V) <p v v - (v u ■ V) V P v ) dxdr 
o Jn 

+v f [ v u ■ (Vn) T • 4>dSdr = [ <0(O) + Po^(0) dx, 
Jo Jr Jn 

for all vector-fields <f> £ C °°([0,T[xn) such that V • <f> = in £1 x [0, T[, and <f> ■ n = on V x [0,T[ 
and scalars ip € C^([0,T[xf!) (In particular the resulting solution satisfies the equations in the 
sense of T>'((0,T) x ft)); the following energy estimate 

Mi^ff +i/ f\\Vv v fdr + u f [v»-(Vn) T -v»dSdT 

\\v u \\ 2 -4- \\n u \\ 2 

< Foil +IIPoll 
2 

is satisfied for all t € [0, T] . 

With this definition we have the following result. 
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Theorem 4.2. Let be given any positive T > ; Vq G L 2 (S7), and Pq G L 2 (f2), suc/i £/ia£ Uq is 
divergence- free and such that Vq ■ n = on V . Then, there exists at least a weak solution of the 
Boussinesq equations © on [0, T\. 

Proof. By simplifying a procedure typical of compressible flows, as in P.-L. Lions |17| and Feireisl 
§ 4], we first consider the following approximate system, with e > 

d t v u ' e - uAv u ' e + {v"' e • V) v u ' e + Vp u ' e = -p u ' e e 3 in fl x (0, T], 

d t p u ' e + (v u ' e ■ V) p u ' e = eAp u ' e in Q x (0, T], 

V-v u ' £ = inOx(0,T], 

v u ' e -n = onrx(0,T], 

cut\v u ' £ xn = on T x (0, T], 

n-Vp"' e = onrx(0,T], 

v u > e (0,x) =v% in ft, 

p"' e (0,a;) =Po in ^- 

By standard techniques one can show that the above system has at least a weak solution 

v u ' e , p u ' e G L°°(0, T; L 2 (ft)) n L 2 (0, T; tf 1 ^)), 
satisfying the energy inequality 

\\v u ' e (t)\\ 2 + ||p^(t) 112 ^ * 



J 

z/' e • (Vn) T • dSdr < 



+ /' / I W' e || 2 dr + e / ||V/9^|| 2 dr 

^112 i ||„i/||2 



r 



2 

As e goes to zero, (and by standard results of compactness) one can find sub-sequences such that 

v v ' e ^v v and p v *^p v in L°°(0, T; L 2 (ft)), 

u"' e in L^Tjtf^Q)), 
p^^p e and v v ' e -> v v in L 2 (0,T;L 2 (ft)), 

for some v u G L°°(0, T; L 2 (S7))nL 2 (0, T; ff 1 (f2)) and p v G L°°(0, T; L 2 (ft)). By standard arguments 
it turns out that (v u , p u ) is a weak solution to ([5]). In particular, being v u ' e strongly convergent and 
p u,e weakly convergent in L 2 (0, T; L 2 (f2)), one has 

f (/' £ • VV') p u ' e dxdr -+ f f (v v ■ VV») P v dxdr, 
o Jn Jo Jn 



and also 



>o Jn 

r-T 



Vp u ' e Vip dxdr -> 0, 



due to the uniform bound for e J Q \ V p u ' e \ 2 dxdT. This implies that the (v v , p u ) is a weak solution 
in the sense of Definition 14.11 and in particular that energy inequality (|17h is satisfied. □ 
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From the existence result we prove now the last result of this paper. 

Theorem 4.3. Let £1 be a bounded smooth open set in R 3 , let v E € H 3 (Q), be a divergence- free 
vector field tangential to the boundary, let also p E € H 3 (Q). Let us suppose that the initial data 
satisfy the following conditions at the boundary: 

(18) ° V 

Vpf(x) = ViGT. 

Let u E , p E G C([0, T]\ H 3 (Q,)) be the unique solution of the Euler Boussinesq equations (|6|), with 
initial datum (v E ,p E ) and defined in some interval [0, T\. 

Let [v v ,p u ) be a weak solution of the Boussinesq equations ([5|) with a divergence free and tan- 
gential to the boundary initial datum Uq £ L 2 (Q) and with pg € L 2 (Q) such that 



K - u E \\ = 0(u) and - p E \\ = Oiy). 



Then, 



sup |K - u E f Loo{0TL2m = 0{u 2 ), sup \\p u - p E \\ 2 L oo {0TL 2 m = 0{u 2 ), 
t£[0,T] y y " te[o,T] y K " 

and 

\\Vu» -Vu E \\ 2 L2 ^ T . L i {n)) =0(u). 

Proof. The proof is based on the same techniques employed before. We multiply © by (v E ,p ) 
and with suitable integration by parts we get 



(19) 



Jo Jn ^ 

+ {v v -V)p u p E )dxdT - [ [ p v p E + v v vf dxdr + u [ [ v v {Vn) T v E dSdr 
' Jo Jn Jo Jr 



v v {t)u"{t) + p v {t)p 1 "{t)dx + 

■t 



vE v o + PoPo dx. 



We then obtain from the equation for (v E ,p E ) 

(20) f [ [vfv v + pf p v + (v E ■ V) // + (v E ■ V) p E v v ) dxdr = [ v E v v Q + p E p v dx. 
Jo Jn Jn 

Next, by multiplying the Euler Boussinesq equations by u E and by the usual integrations by parts 
we get the energy conservation 

(21) \\v E (t)\\ 2 + \\p E (t)\\ 2 = II^H 2 + lk g || 2 
1 ' 2 2 

Then, by adding together (fTT)) - (f2~T]) and subtracting p9|) - (p0|) . we get 

IK*)II 2 + IIpWII 2 



+ u I I Vv v Vvdxdr + / / (v-V)v E vdxdT 
Jo Jn Jo Jn 



2 

I f {v-V)p E pdxdT + v I I v v {VnfvdSdT 
Jo Jn Jo Jr 



< Kl±Mof 
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where 



v := v u — v and p := p v — p E . 



With the same manipulations employed in the previous section we get 



2 rt 
-V 

Jo Jr 



- [ [ (v • V) w E u + (v • V) p E pdxdr - v [ I Av E v dxdr 
Jo Jn Jo Jn 

+ v \ I (curl v E x n) v dS dr. 
Jo Jr 



We then estimate most of the terms as before with 



(v -V)v E v + (v V)p E pdx 



< 



C(||n(t)|| 2 + ||p(t)|| 



v ■ (Vn) T -vdS 



< Cv\\v\\ 2 + -||Vv|| 2 , 



Av E v dx 



<C(\\v\\ 2 + u 2 



where we used that v E G C([0, T]; H 3 (U)). 

To handle the last integral involving curlv^ on F, we observe that the equation for the vorticity 
for the Euler-Boussinesq system implies, along path-lines, 

cuvlv E (X(a, t),t) = V a X(a,t) curl v E (a,0) - [ cuil(p E e 3 )(X(a, a), a) da, 

Jo 

hence now cw:\v E (x) = on T is not enough to have curl-u^ = on the boundary for all positive 
times and a control also on curl^^es) at the boundary is needed. We then observe that p is 
transported by the velocity v E (it solves the equation of continuity) , hence 

p E (X(a,t),t) = p E (a). 

Consequently, we get the following evolution equation for the gradient of p: 

S/p E X(a, t),t) V a X(a, t) = Vp E (a), 

hence 

Vp E (X(a, t),t)= Vp E (a) [V a X(a, t)} '\ 

Since the matrix V a X is non-singular and since path-lines starting at the boundary remain at the 
boundary, it follows that a sufficient condition to have curl(p E es)(x,t) = for all x G T and for 
all t G [0, T] is that of asking (|18p . Under the above assumptions the term Jpfcurlt^ x n)vdS 
vanishes identically, hence we arrive at the inequality 

IK*)H 2 + \\p{t)\? + » f l|V«|| 2 dr < \\u(0)\\ 2 + \\p(0)\\ 2 + C\ f \\u(t)\\ 2 dr + v 2 
Jo 1 Jo ■ 

from which we have the thesis by applying the Gronwall lemma. □ 
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